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I.  INTRODUCTION 

Partial  arc  self-acting  bearing  pads  are  often  used  to  overcome 
rotor-bearing  instability.  In  such  applications,  depending  on  the 
specific  design  requirements,  the  bearing  may  be  constructed  either  of 
a  single  fixed  pad,  of  several  fixed  pads  (multi-lobe  bearing),  or  of 
several  pivoted  or  flexibly  supported  shoes.  The  vital  design  informa¬ 
tion  consists  of  both  the  static  and  time  dependent  fluid  film  forces 
of  each  pad,  which  ultimately  determines  the  load  carrying  capacity 
of  the  bearing  and  the  stability  characteristics  of  the  rotor-bearing 
system. 

An  important  application  of  the  partial-arc  type  bearings  is  the 
self-acting  gas  journal  bearings  operating  in  a  high  ambient  pressure. 
Here  the  lateral  unit  load  on  the  journal  bearings  is  extremely  small 
so  that  Half -Ere quency  Whirl  instability  is  likely  to  occur  for  plain 
cylindrical  journal  bearings.  Under  this  condition,  the  bearing  com¬ 
pressibility  number,  A,  is  always  very  small,  thus,  it  is  convenient 
to  solve  the  lubrication  problem  by  an  expansion  of  the  fluid  film 
pressure  in  terms  of  a  power  series  in  A  [6].  The  first  order  effect 
of  this  expansion  actually  gives  the  same  equation  applicable  to  an 
incompressible  lubricant.  The  higher  order  effects  are  concerned  with 
equations  of  increasing  complexity,  but  are  nevertheless  all  linear. 

Solution  of  these  equations  is  no  simple  matter  because  of  the 
variable  coefficients.  Tanner  ro  has  applied  the  method  of  Galerkin 
to  solve  the  steady  state  Reynolds  equation  for  a  variety  of  bearing 
configurations.  The  present  paper  extends  this  method  to  consider  both 
time -dependence  and  higher  order  effects  of  A.  Under  the  assumptions 
made  herein,  the  results  also  may  be  applied  to  the  calculation  .of 
critical  speed,  rotor  response  [7j,  as  well  as  the  stability  of  flexible 
rotors. 


*Half-Frequency  Whirl  -  A  special  case  of  instability  generally  associated 
with  self-acting  journal  bearings.  This  instability  occurs  when  the  journal 
speed  reaches  a  critical  value.  The  journal  axis  whirls  at  a  frequency  of 
one-half  or  nearly  one-half  of  the  journal  speed  in  the  same  direction  as 
the  journal  rotation.  The  motion  of  the  journal  of  the  journal  axis  can 
be  either  conical  or  cylindrical. 
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II. 


DERIVATION  OF  EQUATIONS  TO  BE  SOLVED 


The  object  of  this  section  is  to  obtain  the  Reynolds  equation 
in  the  form 

L  {U]  =  f’ 

where  L  is  a  linear  operator,  f  is  known  function  and  u  is  the 
desired  solution.  We  may  then  obtain  our  solution  by  Galerkin's 
method.  Thus,  we  consider  the  nondimens ional  Reynolds  equation 

L  {p2}  s  [h3  §t“]  +  a7-  [h3  ar]}{p2 }  =  (hp)  +  2  ar  (hp)]  * 


(±  -f-  ,  z)  =  p2  <0,  ±  L/D)  =  1. 


(1) 


We  make  the  following  three  assumptions: 

(I)  A  is  small  (i.e.,  may  be  used  as  a  perturbation  parameter), 


(II) 


p  does  not  depend  on  time  explicitly;  i.e. , 


<l£  =  <>£  *  +  <*£  •* 

dt  de  be 


+  ^  a  + 

OQ£  oa 


(III)  Time  derivatives  (e.g. ,  €,a,€,a)  Are  small  enough 
so  as  to  make  th^ir  products  negligible. 

Then,  using  Assumption  (I),  we  write 

p(0,z)  =  k2)  Ak  p^k\ 


(2) 


where 


,(o>  ,  1. 


Thus, 


2,„  v  <S  Ak  ^ 


p  (0,z)  =  k£)  A 


n(j)  „(k-J) 

P  P 


(3) 


j=o 


Substitution  of  (2)  and  (3)  in  (1)  gives 

L 


H-ft-  d-3  & +  -k  l*3  4U  p0)  - 


(4) 


from  which,  equating  like  potoers  of  A,  we  obtain 
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jo  4<jyki  -  4R [h3  ^-<pUVk'j>>]+ &  p  ^  *o>»(k'j) 

•  2{sr  +  2-5r}  hp<k  l>i  k  =  1,2. .. .  (5) 

Under  Assumption  (II)  above,  (5)  becomes 


k 

2 

J-0 


fo,p<k-i)}  -  2y- 


21  tR  +  2  2  (e1*1 


r=l 


a 

ae  C'-a 


+  a«  A..,,)]  (hp^ 


aa1 


(6) 


) 


where 

£  W  -  drf /dtr 

k  -  1,2,... 


Vfenote  then,  that  for  k  =  1,  since  p^  =  1,  (6)  becomes 

«■  f(U}  -  R  O-3  £■>  *  ^  ^-)}{p(1)}  -  |  ♦  2  [«  S  *4  1]  . 


(6a) 


For  k>l ,  we  transpose  all  terms  containing  p^\  with  i<k,  to  the  right 
side  of  the  equation  to  obtain 


We  now  represent  each  perturbation  pressure  p  '  as  the  sum  of  its 
steady  state  component  and  its  time -dependent  components  in  the  form: 


.<k>  =  „<k> 


=  P 


00 


k 

+  2 
r=l 


(eWp 


(k) 

rO 


k  -  1,2, 


(8) 


(from  (6a),  we  note  that  for  k  =1: 
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Using  (8)  and  Assumption  (III),  we  may  write 


(j)D(k-j)  „  _<j)  p(k-j)  +  (k-j)  J  ,  M  (j)  +  0Mp(JK 

p  p00  p00  pnn  ^  pr-n  T  u  • 


+  pS>  k*j  (.Wp2“j>  +  «wpS"j)>. 


k  =  2,3. 
j  -  1,2, 

Then  substitution  of  (8)  and  (10)  in  (7)  yields 


+  4  ST  <hp00_1))]  +2t,Vert’h^:l!o+“W<r).l> 

+  a^PQr>)| 

(ID 

for  k  =  2,3, . . . 

We  note  (from  (1))  that  L  is  linear;  i.e.,  L^au  +  bv}“ 

aL  ^u  ^  +  bL  -£v  j-. 

However,  L  acting  on  a  product  is  given  by 

L{uv}=  uL{v}+  vl{u}+  2h3  (f£  g  + 

=  uL-^v^  +  vL^u^+  2M  (u,v) 

(12) 

We  also  now  define  the  ,rknown"  functions  f  (8,z)  by 

TJn(k)l_ 

L1P00  J  foo  ’ 

•  LOc)]..  f(k) 

4"t0  J-  £r0  • 

L  a 

*  J"_(k)”L _  r (k) 

L  jP0r  J  ~  f0r  ’ 

(13) 

5 


noting  that  any  ^  occuring  in  an  expression  for  f^k^  has  J<k. 
A.  Steady  State  Equations 

Consideration  of  (11),  and  use  of  the  relations  (12)  and 
(13)  gives 


fOO  _  &  ,hn(k-lK  . 
foo  -  W  (hpoo  > 


k-1 

E 

j=l 


f(k-J) 

00 


+  P 


(k-J) 

00 


2M  (p 


(j) 
00  ’ 


p00  ' 


> 


(14) 


for  k  =  2,3, .. . 


B.  Time-Dependent  Equations  (k  =  2,3,...) 

Again  considering  (11),  matching  like  e-  and  Of-  derivative 
terms,  and  using  the  relations  (12)  and  (13)  gives,  after  some 
manipulation: 


r(  „(k)\  _  f(k)  _  *  (k-1).  ,  /.  (k-1  )v 

Ll  pio  /  “  fio  -  2  &  (hp  00  >  +  5e  (hpio  ) 


_(r)  f(k-r)  (k-r)  (r)  (r)  (k-r). 

‘  Z,  P00  f  10  +  P10  f00  +  2M  (P00  *  p10  > 
r=l 


L(  p(k)l  =  f(k)=  2hp(k‘1>  +  ^  -  Z  (viT)  f(k'r)+  0(k'r)  f(r) 

lPj0;  rj0  Zhpj-l,0  +  ^0  <hpj,0  >  r“x  (p00  fj0  +pjo  f00 

J  -  1.2,3 . k-1  (15) 


+  »  cp£>.  . 


w*-l) 


k-1 ,0  ’ 


Tfn(k)L  f<k>  -  9  *  /i__ (k— 1 ) \  .  *  r.  (k-1).  k:X  (r)  .(k-r) 

Llpoi  r  foi  -  2  55  (hp  00  >  +  $e  (hpoi  >  ’  ^  (poo  f  01 


.  (k-r)  (r)  ?M,  (r)  (k-r). 

+  P01  f00  +  2M(poo  ’  P01  >  ’ 

k-j 


L(p(k)l  -  f(k)-  2hD(k_1)  +  fhD(k_1)*  V  rn(r)  f(k'r>  +  n(k'r)f(r) 

Llp0j  f  ~  f0j  -  2hp0,j-i  +  *e  (hp0,j  >  rfx  (p00  f0j  +  p0j  foo 


+  2«<4o)’  Cr)>  • 


j  *  2,3,... ,k-l. 


,(k)\  =  f(k) 
0k 


{C}- 


uPr'j.=  fxr'  -  2hP 


(k-i) 

0,k-l. 
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C.  Static  Stability  Derivatives 

We  are  also  interested  in  obtaining  equations  for  the 
derivatives  of  the  steady  state  pressures  with  respect  to 
e  and  a.  Therefore,  we  consider 


We  assume  that: 


(1)  All  necessary  derivatives  exist,  and 

(2)  h  is  independent  of  z. 


Then,  differentiating  the  above  expression  by  e,  we  obtain 


dp, 


00  1 
oo 

e 


=  (k) 


*f(k> 

00 

1e 


-3  (' 


I  f(k>  ,  ah 


dp, 


(k) 

00 


-3h 


2  Sp00 


h  00 
(k) 


a2h 


o6  Oe 

Similarly,  a  differentiation  by  a  yields 


L  < 


<k) 


Sp00 

15“ 


:  (k) 
■  a 


la 


-3h 


2  Sp00 


(k) 

OC 

W 


d2h 

aaae 


ae 


VV  V 


-3  (—  f^  -h  ^ 

'■h  oo  5e 


dp, 


(k)  . 

00  %  oh 

0  '  aa 


Thus,  for  any  power  k  of  the  perturbation  series  (2)  for  p, 
we  have  obtained  a  set  of  partial  differential  equations  for 
the  steady  state  pressure,  the  time -dependent  pressure  com¬ 
ponents,  and  the  derivatives  of  the  steady-state  pressure 
with  respect  to  e  and  a  all  of  which  are  in  the  form 


L{u}=  f 

suitable  for  numerical  solution  by  the  Galerkin  method. 
We  summarize  these  below  for  k  ■  1  and  2. 


k  -  1: 


(16) 


(17) 


{.si  ■ 
{'«)• 


and 


Is?)- 

( 


E(i) 

C10 

E(1) 

C01 


p(D 


2  Se 


2  ir  * 


(9) 


(l-3h 


,(D  2 

2  dp00  ,  b\ 


3  ah  dh 


)U  u  UU  uu  mU 

Sed9  h  d9  de 


( 1 ) 

2  dp00 


bQ 


(1) 

00 


(1) 


n-fi-J-  a 

k  =  2: 


(1) 


(l-3h 


2  dpQQ  .  b\  3  dh  ah  ,2  ^p00 


"a<r  'a^ae '  h  ae  aa 


ae 


(2)  (2)  •  (2)  ..  (2)  &  (2)  a  (2) 
p00  +  €  pm  +  6  Pon  +  a  Pm  +  a  P02  * 


'10 


20 


'01 


'02 


), 

(18) 


-).  (19) 


where 


T  r  (2)1  ad  _  _a_  a).  (i)  (i)  _  ,  (i)  (Di 

L  lP00  J  =  foo  ~  ao  <hpoo  '  poo  foo  M(poo  ,P00  J 


*„(!) 

,  dp00  ,3 

h  "ae"  •  h 


i ■ 


{plo*}  -  «{?  ■  §t  +  2  **-*»£'>  -  f 


,(1). 

10 


_a_ 

ae 


/l), 

00 


(l)f(l) 
00  £10 


(!),(!)  2Mfo(l)  d(Dn 

■p10  foo  2M(p00  ,pio  '• 


:(2)  _ 
L20 


da  =  2hP 


(1) 

10 


f  (2)1  ,(2)  a  (Ik  .  -  d  ,.(lk  . 

|poi  j  ~  foi  =  ae  (hpoi  ^  +  25a  (hpoo  ^  poo  foi 


ADAD  _  2M  ,_(U  AD, 
"p01  foo  2M  'P00  ,P01 


{’$} 


f(2) 

02 


2hp 


(1) 
01  ’ 


Y  (20) 


(21) 

(22) 


(23) 


(24) 


8 


and 


(25) 


(26) 
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III.  SOLUTION  BY  GALERKIN'S  METHOD 

In  the  present  report,  we  outline  briefly  the  steps  followed 
in  solving  our  equations  by  Galerkin's  method.  This  method  is 
equivalent  to  the  Rayligy-Ritz  procedure  in  any  problem  where  a 
minimum  principle  applies,  as  in  this  case.  A  complete  discussion 
is  given  in  [2]. 

We  consider  the  partial  differential  equation 

1  {-}  ■  -Is <"3 1  >  +  -fr  <»3  tr  >  ■  £  «•*>•  <”> 

We  assume  that  (See  Figure  1.) 

h  *  1-e  cos  (0-a).  (28) 

As  boundary  conditions,  we  have 

u  (±  |  ,  z)  =  u  (0,  +  |)  =0.  (29) 

From  the  above  it  follows  that  u  is  an  even  function  of  z,  i.e.. 


f  (0,z)  =  f  (0,-  z) 


(30) 


We  assume  an  approximate  solution  u  (0,z)  of  the  form 


N  M  _ 

u  (0,z)  -  Z  Z  (z)  2  c  ®  (0), 


n*l 


m»l 


mn  m 


(31) 


where  the  functions  Zr  (z)  and  (0)  are  1  inear ly,  inde¬ 

pendent  elements  of  a  complete  set  (i.e.,  a  set  complete  in  the  sense 
of  the  Weierstrass  approximation  theorem  [3]  ,  and  satisfy  the 


boundary  conditions  (29).  In  our  case,  we  choose 
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Zn(z)  -  cos  pSill  z]  , 

(32) 

®m  (9)  "  8ln  [  f  6  +  f]  • 


The  c  are  undetermined  coefficients, 
mn 

Then,  following  the  procedure  outlined  in  [2]  ,  we  determine  the 

coefficients  c  from 
mn 


P/2  L/D  - 

2  4/2  4  LL  Vz)  V0)dz  de  -  0, 


(33a) 


for  V  =  1,2,. ..,N;  p  =  1,2, ...,M.  Or, 

z  z  ;P/  / /  [z  df" (h3  +  h3®  — *y-\  zv©  dz 

a-1  m-1  -P/z  0  L  n  d®  d9  m  dz2  J  V  u 


d2Z 


d6 


p/2  _  L/D 

-  0  I 


f(6,z)Zv  dz  d0, 


W  p  U 

for  V  =  1,2,. ..,N;  p  =  1,2,. . . ,M. 


(33b) 


Ordinarily,  to  obtain  the  coefficients  c  ,  we  would  have  to 

mn 

invert  a  matrix  of  order  MN  x  MN.  However,  with  our  choice 
of  Z^(z)  and  the  assumption  on  h  (Eq.  28),  the  orthogonality 
of  the  cosine  functions  requires  that  n  =  V,  and  (33b)  becomes 
partially  uncoupled.  Therefore,  we  have  instead  N  matrices 
each  of  order  M  x  M  to  invert. 


In  the  preceding  section,  we  derived  the  forcing  terms  f 
corresponding  to  each  pressure  component  and  the  steady  state 
pressure  derivatives  with  respect  to  e  and  a.  We  now  corres¬ 
pond  to  each  of  these  a  set  of  coefficient  matrices: 


f(k)^ 

rs 


€ 


,00  = 


'  rs 


=  Z 
n=l 


M 

Z 

m=l 


dp 


(k) 

00 


de 


N 
=  Z 
n=l 


M 

Z 

m«l 


(34a) 


(34b) 
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(k) 


a 


da 


N  M  ,  > 

Z  Z  ;  9  (0)  Z  (z). 


i  i  mn 

n=l  m=l  (x 


n 


(34c) 


The  uncoupling  in  z  described  above  permits  us  to  write  the 
systems  to  be  inverted  as  follows,  for  each  value  of  n  (n  «1, 


M  (k) 

Z  cl  '  A 
,  mn  pm;n 

m«l  rs  ^ 


(k)  =  _|£D  ^  g|  j  1  f(k)  z  dz  dQ 

|jn  pL  *p/2  p  ■'0  rs  n 


M 


Z  c<k>  A 

mn  pm;n 


m*l 


p  =  1,2 . M, 


,<*> 


2flD  rP/2 

PL  ip/2 


(35a) 


kn  Vo 


L/D 


f(k)  Z  dz  d0; 
e  n 


P  =  1 , 2 , . . . ,M, 

2  c(k)  A  -  .W-IP  fn  6  /L/V>  Z  dzde; 

m-1  “  imin  an  P  -P/2  4  0  €  n 

p  =  1,2, ...M. 

In  each  of  the  above,  the  elements  of  the  matrices  A  are 

pm  ,n 

identical,  and  are  given  by 

3 


pm;n 


<v +  x»> +  \  jfc  { aj  c°s 


1  - 


2(k 2  +  X2)1 
m  n 


j' 


x 


,ln  (.18/2)11  _  . 

#)  -  4„2  JJ 


■m {i 


2  fa.  sin  ja  fn2  +  k2  +  2\2  -j2"|  [ - 2  -C-S-^ 

Lm  4  n  J[|^-mV 


(35b) 


(35c) 


LM21 _ 

2]  2-(2mp)2J 


1, 


m  t  p 
m  +  p  odd 


(36) 


-mp 


jfL  {*j ce>* Ja  +  V +  2xn--)2]|^H?^r 


(2mp) 


m  t  p 
m  +  p  even 
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where 


*  ^  3  2 . 

—  (1  +  —  e  )  , 

.  M  2 


..a  jl.  3 

2  a  € 


K  =  mx/f3  , 


\ 

n 


,  2n-l.  flD 
{  2  }  L 


Note  that  the  A  matrices  are  symmetric  in  p  and  m.  They 
also  contain  some  removable  singularities.  The  limiting  values 
of  terms  in  (36)  are  given  as  follows: 


lim 

x— ►  2m 

lim 

x-*|p  +  m| 


lim 

x-*-||i  +  m| 


■ 


sin  Qtx/2)  1  =  (- )mrt  t 


2  /  2 
x  -  4m 


j 


;  n  =»  m 


cos  (itx/2)  1 

-  *  st"[f  1“  ±  ”|]  fu  >»  1 

2  2  2.2  ..  . 2  [  ~ 
[x  -m  -p  )  -(2m p)  J 

16  pm  |p  +  m|  *  1  p  +  m  oddT 

[sin  (Jtx/2)  { 

«co,[fp  +  .j  fj.>»  1 

}(x2-m2-(i2)2  -(2mp)2J 

-  ±  16  pm  |p  +  m|  ^+”evenf 

> 

>(37) 


The  computational  procedure  to  obtain  pressure  is  then  as 
follows : 

1.  Calculate  the  A  matrices  from  (36)  and  invert  them. 

pm;n 

2.  For  the  appropriate  pressure  component  obtain  the  force- 
ing  term  from  Eqs.  (9),  (14)- (17).  These  are  given  in 
more  detail  for  k  =1,2  in  Eqs.  (9),  (18)- (26)  on  pages 


3.  Obtain  the  matrix  of  4>n  vectors  from  Eqs.  (35). 

4.  Obtain  the  coefficients  c  by  multiplying  the  ®  vectors 

mn  n 

by  the  corresponding  inverses  of  the  A  matrices. 

J  r  |im;n 

5.  The  pressure  component  is  then  found  using  Eqs.  (34),  (8) 
and  (2). 
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IV.  CALCULATION  OF  FORCES 


The  dimensionless  forces  along  and  perpendicular  to  the  line 
of  centers  are  denoted  FR  and  FT  respectively,  and  are  given  by 

P/2  L/D 

F_  ■  2  /  cos  (0-OC)  /  p  dz  d0, 

R  -P/2  0 

P/2  L/D 

F_  *  -  2  /  sin  (0-Q!)  /  p  dz  d0, 

T  -p/2  0 


Where  p  is  the  pressure  component  or  its  derivative  with  respect 
to  e  or  Ct.  Substituting  the  appropriate  summation 


P 


N  M 

Z  Z  c  ©  Z 
mn  m  n 

n=l  m=l 


in  (38),  we  may  then  integrate  directly  to  obtain 


8 

L 

N 

v 

<->n 

M 

y 

c  r  -  1 

n 

D  *  ‘ 

Li 

n=l 

• 

2n-l 

Li 

m*l 

“  L»2-<p/»)2J 

-cos  ^  cos  a,  (m  odd) 


X 


sin  sin  OC,  (m  even) 


cos  ( p/2) sin  a, 
sin  (p/2.)cos  a, 


(38a) 


(38b) 


(39a) 


(39b) 


V. 


EXTRAPOLATION  PROCEDURE 
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By  use  of  Richardson's  extrapolation  procedure  (see  e.g.  [4]  ) ,  and 

Tanner's  use  *  of  extrapolation  in  [1]  ,  we  may  use  results  computed 

with  different  values  of  M  and/or  N  to  obtain  a  closer  approximation 

to  the  true  solution.  We  will  assume  that  the  z  dependence  of  our 

solution  can  be  obtained  with  very  few  functions  Z  .  Henceforth. 

n 

we  consider  N  fixed.  Then,  if  we  use  two  values  of  M,  we  obtain  the 
extrapolation  formula 

F  -  (M*  E2  -  M*  Ex)  /  (M2  -  mJ) ,  (40) 

where 

F  *  extrapolated  result, 

*  computed  result  for  M  *>  M^;  i  ■  1,2. 


It  is  of  interest  to  include  a  short  discussion  of  error  growth 
in  the  calculation  of  f^  for  successive  values  of  k.  We  assume 
truncation  error  per  calculation  fixed.  (That  is,  we  assume  M  and 
N  are  fixed.  Obviously,  the  greater  M  and  N  to  begin  with,  the  lower 
the  truncation  error,  from  the  Weierstrass  theorem).  Thus,  for  a 
given  k,  we  have 


P 


00 


+  € 


(k) 


9 


where 

p^»  calculated  pressure  at  the  k~  stage  of  computation, 

(k)  .  .th  .  . 

p  *  true  pressure  at  the  stage, 

(k)  .  th  ^ 

€  *  error  at  the  stage. 


(k  +  1)  (k) 

To  calculate  fv  '  requires  derivatives  of  p  .  But,  in 

M  ,  it  is  shown  that  the  errors  in  the  first  derivatives  of  the 

solution  function  are  of  the  same  order  as  those  in  the  solution 


In  [lj  ,  Tanner  used  a  three-term  extrapolation,  wherein  the  first  term 
obviously  was  not  good.  If  he  had  used  his  two  better  approximate  values 
in  an  extrapolation  of  the  form  (40)  above,  his  results  would  have  been 
much  closer  to  those  he  used  as  a  standard. 


function  itself.  Thus,  it  can  be  seen  that  while  error  does  increase 

with  k,  the  process  is  essentially  stable  in  that  all  errors  in  the 

(kl 

(k  4-  1)  stage  are  of  order  e'  ,  and  therefore,  higher  values  of  k 
can  be  computed  if  the  initial  truncation  errors  are  kept  low. 


16 


VI.  SAMPLE  RESULTS 

A  digital  computer  program  has  been  written  to  perform  the 
calculations  indicated  above  for  k  =  1.  Currently  being  checked 
out  is  a  subroutine  to  extend  these  results  to  the  case  k  =  2. 

Some  computer  program  results  are  given  herein.  The  physical 
inputs  were  as  follows: 

Shoe  angle  (3  ■  120° 

Angle  of  line  of  centers  Q!  =  60° 

L/D  -  1.0 

Eccentricity  ration  e  *  0.01,  0.8 

Results  were  computed  with  N  ■  2,  and  with  two  values  of  M:  4  and  6. 
A  tabulation  follows: 


Force 

€  ■ 

0.01 

i _ t-°-8 

M 

M 

Component 

4 

6 

OO 

(Extrapo¬ 

lation) 

6 

00 

(Extrapo¬ 

lation) 

F  ^  m 

Roo 

.lFd) 

2  R01 

0.00187 

0.00187 

0.00187 

1.367 

1.490 

1.589 

r™  m  ' 

Aoo 

1  p  (1) 

2  V 

0.00384 

0.00384 

0.00384 

1.111 

1.140 

1.164 

fr(1) 

R10 

0.3453 

0.3467 

0.3479 

5.533 

6.865 

7.931 

,U) 

Aio 

0.3732 

0.3730 

0.3728 

3.417 

3.726 

3.973 

0.1899 

0.1898 

0.1897 

9.397 

11.081 

12.428 

oo 

1-  P*1* 
^  Too 

0.3869 

0.3874 

0.3878 

5.569 

5.859 

6.091 

h  fr1} 

^  Roo 

0.00169 

0.00170 

0.00171 

-2.030 

-1.693 

-1.423 

S  Too 

0.00178 

0.00177 

0.00176 

-1.861 

-1.807 

-1.764 
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A  plot  of  the  steady  state  forces  and  pressure  profile  is 
given  as  Figure  1.  The  above  results  have  been  compared  with 
those  obtained  in  p] ,  and  good  agreement  was  found. 


CONCLUSIONS 


1.  The  approach  described  herein  provides  a  fast,  inexpensive, 
accurate  method  of  obtaining  partial  arc  and  pivoted  shoe 
bearing  data. 

2.  It  operates  over  a  wide  range  of  conditions  (including  high 
values  of  e),  over  which  difficulty  is  encountered  in  other 
(e.g. ,  iterative)  approaches. 

3.  The  expressions  given  in  the  present  report  include  compress 
sibility  effects  of  all  orders  (i.e. ,  powers  of  A).  The  computer 
program  can  be  extended  to  obtain  these  higher  order  effects. 


RECOMMENDATIONS 


A.  Computer  Program 

It  is  recommended  that: 

1.  The  computer  program  be  extended  to  include  second-order 
effects  (k  =  2). 

2.  The  program  be  modified  so  that  Richardson's  extrapolation 
procedure  is  performed  automatically  by  the  computer. 

3.  Programming  be  included  so  as  to  perform  numerical  quadratures 
for  higher  compressibility  effects  (values  of  k),  since  analyti¬ 
cal  procedures  are  not  practical. 

B.  Applications 

1.  The  computer  program  should  be  used  to  generate  design  charts 
for  partial  arc  gas  bearings  at  high  ambients  to  check  and  sup¬ 
plement  presently  available  data. 

2.  The  dynamic  force  derivatives  calculated  by  this  computer  pro¬ 
gram  can  be  incorporated  into  rotor  dynamics  analysis  to  study 
critical  speeds  and  frequency  response. 

3.  The  static  and  dynamic  force  derivatives  can  be  used  to  calcu¬ 
late  the  threshold  speed  of  instability  of  rotors  supported 

in  partial  arc  bearings. 

4.  This  computer  program  can  be  used  as  a  subroutine  in  a  larger 
program  to  study  static  and  dynamic  characteristics  of  multi- 
lobe  and  pivoted  shoe  bearings. 
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Lynn,  Maaeechueette  l 


Reaearch  Laboratories 
Ganaral  Motors  Corporation 
General  Motora  Technical  Cantar 
12  Mila  and  Mound  Roada 
Warren,  Michigan 
Attn:  Mr.  E.  Roland  Makl, 

Mechanical  Development  Dept.  1 

A.  C.  Spark  Plug  Divleion 
General  Motora  Corporation 
Mllwaukea  1,  Wit cons in 

Attn:  Allen  Knudsan  1 

Mr.  Walter  Cerow 
Kearfott  Diviaion 
Ganaral  Precision  Inc. 

1150  McBride  Avenue 

Little  Falla,  New  Jersey  1 

Grumman  Aircraft  Engineering  Corp. 
Bethpage,  Long  Ialand,  New  York 
Attn:  Mr.  David  W.  Craig,  Jr. 

Mechanical  Design  Section 


Engineering  Department  1 

Hydronautica ,  Incorporated 
200  Monroe  Street 

Rockville,  Maryland  1 

International  Business  Machine  Corp. 
Research  Laboratory 
San  Jose,  California 

Attn:  Dr.  W.  E.  Langlois  2 

Mr.  B.  A.  Napier,  Engr.  Dlv.  Mgr. 

Lear,  Incorporated 
110  Ionia  Avenue  NW 

Grand  Rapids,  Michigan  1 

Mr.  L.  R.  Barr,  V.  P. 

Lear-Romec  Division 
Abbe  Road 

Elyria,  Ohio  1 


Dr.  Calus  G.  GosCzel,  D/53-30 
Bldg.  201,  Plant  2,  Palo  Alto 
Lockheed  Missiles  4  Space  Co. 

P.  O.  Box  504 

Sunnyvale,  California  1 


Dr.  J.  S.  Auamaa 
Litton  Syatams,  Inc. 

5500  Canoga  Avanua 

Woodland  Kills,  California  l 

Mr,  Don  Moora 
Litton  Systems 
3500  Canoga  Avanua 
Woodland  Hills,  California 

Mr.  A.  N.  Thomas 
Astro  Divleion 
Tha  Marquardt  Corporation 
16355  Saticoy  Street 

Van  Nuys,  California  1 

Mr,  Kendall  Perkins 
Vice  President  (Engr) 

McDonnell  Aircraft  Corporation 
Lambert  St.  -  St.  Louis 
Municipal  Airport 
Box  516 

St.  Louis  3,  Missouri  1 

Dr.  Beno  Sternllcht 

Mechanical  Technology  Incorporated 

968  Albany-Shaker  Road 

Latham,  New  York  3 

Mr.  J.  W.  Lower 

Chief,  Engineer-Inertial  Covenant* 

Honeywell  Aero  Diviaion 
2600  Rldgway  Road 

Minneapolis,  Minnesota  l 

Mr.  Carl  F.  Graeasar,  Jr. 

Director  of  Reaearch 

New  Hampshire  Ball  Bearings,  Inc. 

Peterborough,  New  Hampshire  1 

Mrs.  Alice  Ward,  Librarian 

Norden  Division  of  United  Aircraft  Corp. 

Helen  Street 

Norwalk,  Connecticut  1 

Autonetlca  Division  (Librarian) 

North  American  Aviation  Inc. 

9150  East  Imperial  Highway 
Downey,  California 

Northrop  Corporation 
Noralr  Division 
1001  East  Broadway 
Hawthorne,  California 

Attn:  Technical  Information,  3123  1 


Mr.  N.  L.  Sinmons 
Nortronics 

A  Diviaion  of  Northrop  Corp. 

500  Eaat  Orange thorpe  Avenue 

Anaheim,  California  1 

Nortronics 

A  Division  of  Northrop  Corporation 

100  Moras  Strsat 

Norwood,  Massachusetts 

Attn:  Mr.  E.  L.  Svalnaon,  Tech.  Aaat. 

Precision  Products  Department  1 

Pratt  6  Whitney  Aircraft 

Division  of  UAC  -  CAN EL 

P.  0.  Box  611 

Middletown,  Connecticut 

Attn:  Librarlen  1 

Library,  Bldg.  10-2-5 

Radio  Corporation  of  America 

Camden  2,  New  Jersey  1 

Mr.  Robert  S.  Slegler 
Rocks tdyne 

Nucleonics  Subdivision 
6633  Canoga  Avenue 

Canoga  Park,  California  2 

Ryan  Aeronautical  Company 
Attn:  Engineering  Library 

Lindbergh  Field 

San  Diego  12,  Celifornia  1 

Jack  &  Heinte,  A  Division  of 

Tha  Slegler  Corporation 

1725  Eye  Street,  Suite  505 

Washington  6,  D.  C.  I 

Mr.  Paul  A.  Pitt,  Vice  President 
Engineering  &  Research 
Solar  Aircraft  Company 
2200  Pacific  Highway 

San  Diego  12,  California  1 

Mr.  W.  G.  Wing 

Sperry  Gyroscope  Company 

C-2 

Great  Neck,  New  York  1 

E.  A.  Bablay 

Union  Carbide  Nuclear  Company 
Poet  Office  Box  ? 

Oak  Ridga,  T anna • sea  1 
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Mrs.  Florence  Turnbull 
Engineering  Librarian 
Sperry  Gyroscope  Company 

Great  Neck,  New  York  1 

Sunda trend  Aviation-Denver 
2480  West  70th  Avenue 

Denver  21,  Colorado  1 

Mr.  James  W.  Salasai,  President 
Turbocraft,  Inc. 

1946  S.  Myrtle  Avenue 

Monrovia,  California  1 

Universal  Match  Corporation 
Avionics  Department  Technical  Library 
4407  Cook  Avenue 

St.  Louis  13,  Miaaourl  1 

Waukesha  Bearings  Corporation 
P.  0.  Box  346 
Waukesha,  Wisconsin 

Attn:  Mr.  J.  M.  Gruber,  Ch.  Engr.  1 

Mr.  John  Boyd 

Westlnghouse  Electric  Corporation 
Research  Laboratories 

East  Pittsburg,  Pennsylvania  1 

Mr.  H.  Walter 
Director  of  Reaearch 
Worthington  Corporation 

Harrison,  New  Jersey  1 

Dr.  W.  A.  Gross 
Ampex  Corporation 
934  Charter  Blvd. 

Redwood  City,  California  2 

Professor  J.  Modray 
Department  of  Mechanical  Engineering 
Union  Collage 
gchanactady  8,  Haw  York 


1 


